Inward and Outward Node Accessibility in Complex Networks as Revealed by 

Non-Linear Dynamics 
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In this work, the outward and inward accessibihties of individual nodes are defined and their po- 
tential for application is illustrated with respect to the investigation of 6 different types of networks. 
The outward accessibility quantifies the potential of an individual node for accessing other nodes 
through random walks. The inward accessibility measures the potential of a given node of being 
accessed by other nodes. Both the inward and outward accessibilities are measured with respect to 
successive time steps along the walks, providing an interesting means for the characterization of the 
transient non-linear dynamics of accessibility. Self-avoiding walks are considered here because they 
are more purposive and necessarily finite (unlike traditional random walks). The results include the 
identification of the fact that the inward values tend to be much smaller than the outward values, the 
tendency of the inward accessibility to be highly correlated with the node degree while the outward 
values are mostly uncorrelated with that measurements, the distinct behavior of the accessibility 
in geographical networks, the dominance of hubs in scale free networks, as well as the enhanced 
uniformity of the accessibilities for the path-regular model. Also investigated was the possibility to 
predict the accessibility of a given node in terms of its respective degree. The concepts of inward and 
outward accessibility, as well as the several obtained results, have several implications and potential 
for applications to several real-time problems including disease spreading, WWW surfing, protein 
interaction, cortical networks and network attacks, among others. 

PACS numbers: 05.10.-a, 05.40Fb, 89.70.Hj, 89.75.k 



'A map of Esmeralda should include, marked in differ- 
ent colored inks, all these routes, solid and liquid, evident 
and hidden. ' (Invisible Cities, I. Calvino) 



I. INTRODUCTION 

One of the main current issues in complex systems, 
defined by an intersection between dynamic systems and 
complex networks, regards the interaction between struc- 
ture and dynamics (e.g. [JQ)- A great deal of the cur- 
rent attention has been concentrated on linear dynamics, 
especially synchronization (e.g. 0), allowing interesting 
results on how the several structural features of networks 
affect this type of dynamics. Yet, while linear dynamics 
involving diffusion (related to traditional random walks) 
and synchronization are of extreme importance to a wide 
range of important natural and artificial systems, while 
also providing a first order approximation to several non- 
linear dynamics, growing attention is poised to focus the 
extension of the structure-dynamics investigations to in- 
trinsically non-linear dynamical systems. Previous re- 
lated approaches include the studies of associative mem- 
ory capacity in Hopfield networks with connectivity is 
underlain by complex networks [1, |3| , the study of Ising 
and Potts dynamics in complex networks (e.g. @, @, 
opinion formation (e.g. iSj, j3j), among others. 

Because of their simplicity of implementation in com- 
puter systems, random walks stand out as a particularly 
interesting approach to investigate the relationship be- 
tween structure and dynamics in non-linear dynamical 
systems. While many works have addressed traditional 



random walks in complex networks (e.g. [13, El, [13, [H, 
[13, [ill), relatively fewer articles have reported studies 
related on non-linear random walks (e.g. [H, [13, [3, [H, 
[Tgl . [20|). More recently, the transient evolution of self- 
avoiding walks on several types of networks was quanti- 
fied in terms of diversity entropy signatures [2l| . 

The current work extends and complements the pre- 
vious investigations reported in |21|] by considering the 
accessibility of each individual node in the network as 
quantified by the frequencies of visits along self-avoiding 
walks initiating at all the other nodes and introducing 
the concept of inward accessibility. While the concept 
of diversity adopted in that work related to the diversity 
of nodes (not edges) reached by random walks, here we 
use the term accessibility in order to quantify the esti- 
mated number of accesses received or made by specific 
nodes along successive steps of the walks. Two situa- 
tions are considered regarding node accessibility: (i) how 
many nodes can be potentially accessed after starting a 
self-avoiding random walk from a specific node; and (ii) 
how many accesses a particular node can potentially re- 
ceive from self-avoiding random walks starting from all 
other nodes. These two properties are henceforth called 
outward and inward accessibilities. In both cases, it is 
interesting to quantify objectively the accessibility with 
respect to each subsequent number of steps. The current 
article proposes intuitive and sound means to quantify 
the inward and outward accessibilities. The definition 
of these two features naturally motivates the investiga- 
tion of possible asymmetries of accesses, e.g. quantified 
by the ratio between the outward and the inward acces- 
sibilities. We also consider networks from 6 families of 
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models presenting distinctive structures, namely Erdos- 
Renyi, Barabasi- Albert, Watts-Strogatz, a simple geo- 
graphical type of network, as well as two knitted net- 
works (PN and PA) introduced recently [23]. Several in- 
teresting results are reported and discussed regarding the 
inward and outward accessibilities, which are shown to 
present properties intrinsic to each type of network. An- 
other main result is that the outward accessibility tends 
to be substantially larger than the inward accessibility. 
In order to try to relate the accessibilities (a dynamical 
property of the networks) with structural features, we 
consider the possible correlations between the accessibil- 
ities and the degree of each node in the network. It is 
shown that though a strong positive correlation is verified 
between the inward accessibility and degree (especially at 
the later steps of the walks), while the outward accessi- 
bility can only be predicted from the node degree (i.e. 
strong correlation between these two measurement) for 
the very first step of the self-avoiding walks. 

The article starts by presenting the basic concepts, in- 
cluding the motivation and definition of the inward and 
outward accessibilities, and proceeds by presenting and 
discussing the results with respect to individual node ac- 
cessibility and correlations between the latter and the 
node degree. 



II. BASIC CONCEPTS 

This section presents the basic concepts and methods 
used in this article regarding complex network represen- 
tation and characterization, complex network models, as 
well as the definitions of inward and outward accessibil- 
ities and a simple computational method for their esti- 
mation. 



A. Network Basics 

Complex networks can be understood as graphs with 
particularly intricate structure. A unweighted, undi- 
rected complex network can be completely represented 
in terms of its adjacency matrix K of dimension N x N, 
so that each edge extending from node i to node j implies 
K{i,j) = K{j,i) = 1, with K{i,j) = K{j,i) = being 
otherwise imposed. The immediate neighbors of a node 
i are those nodes which are linked to i through a single 
respective edge. Two nodes are said to be adjacent if one 
is the immediate neighbor of the other, and vice versa. 
Two edges are said to be adjacent whenever they have at 
least one of their extremities attached to a same node. 
The degree of a node i corresponds to the number of its 
immediate neighbors. Nodes with degree one are called 
extremity nodes. 

A walk is a sequence of adjacent links, starting at an 
initial node and proceeding along successive steps h. A 
self-avoiding walk is a walk which never visits any node 
or edge more than once. A path is the subgraph asso- 



ciated to a self-avoiding walk. Because of their intrin- 
sic nature, paths represent the most economic way, i.e. 
involving the smallest number of edges, to interconnect 
a set of nodes. In addition, unlike traditional random 
walks, self-avoiding walks necessarily terminate (i.e. the 
moving agent always reaches a point from which it can 
no longer proceed because of lack of unvisited edges or 
nodes) . The nodes where the self-avoiding walks end are 
henceforth called termination nodes. In the present work, 
the moving agents are assumed to stay at the terminal 
nodes after reaching them. This choice allows the con- 
servation of the number of moving agents along all time 
steps. See [2l[ for an additional discussion of such an 
assumption. 



B. Complex Network Models 

Six different types of networks [J, [3, [H, [H, [13, ^ 
are considered in this article, exhibiting markedly dis- 
tinct structural properties. The Erdos-Rcnyi (ER) (see 
also [1^) are obtained by undirectionally connecting N 
initially isolated nodes with a constant probability 7. 
The average degree of this network is (k) = {N — 1)7. A 
Barabasi- Albert (BA) network can be obtained by start- 
ing with mO nodes and incorporating additional nodes, 
one at a time, with m edges, which are attached to the 
previous nodes with probability proportional to their re- 
spective degrees. The average degree of this network is 
(k) = 2m. The Watts-Strogatz networks (WS) can be 
generated by starting with a completely regular network, 
henceforth a linear regular structure, where each node is 
connected to its neighbors at each side. Afterwards, a 
percentage of edges (in our case lOof existing edges) are 
randomly rewired. The geographical type of network con- 
sidered in this article is obtained by distributing N nodes 
uniformly along a two-dimensional space and connecting 
every par of nodes with distance is smaller or equal to a 
maximum distance dmax- The two knitted networks (2^ 
used in this article are the path-transformed BA network 
(PA) and the path-regular network (PN). The former is 
obtained by transforming the star-connections in a BA 
network into respective paths [2^, [23| , so that the result- 
ing network contains a power law distribution of path 
lengths. The latter type, i.e. path-regular network, is 
particularly simple and is obtained by applying the fol- 
lowing procedure several times: (a) start with N isolated 
nodes; (b) choose an initial node and make it the current 
node; (c) select randomly, but without repetition, a new 
node among the remaining nodes and connect it to the 
current node; (d) repeat step (c) until all nodes are con- 
nected through a self-avoiding path. This kind of network 
has been found (2^ . [28j to exhibit remarkable uniformly 
with respect to a large series of measurements which is 
much higher than that presented by the ER model. 

All networks are grown so as to have approximately 
the same average degree {k) and number of nodes N . 
We assume K{i,i) = in all networks considered in this 
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work. Also, only the largest connected component of each 
network is considered. However, because of the relatively 
large the average node degree adopted in this work ((fc) = 
2m — 6), most nodes will be typically contained in the 
principal connected component. 



C. Inward and Outward Accessibility 

Consider the simple but particularly important net- 
work shown in Figure [TJa), involving a central node 1 
connected to 7 surrounding nodes, totaling N = 8 nodes. 
Let a moving agent start a non-preferential self-avoiding 
walk from node 1. After one step, the agent may be at 
any of the surrounding nodes 2 to 8. If we repeat this 
walk M times and count the number of times Q{2) that 
the agent reaches node 2, we can define the rate of ac- 
cesses to that node as 



(a) 



(b) 








A(l,2) 
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More specifically, in the case of Figure [l] the above 
limit converges to Pi (1,2) = 1/7. Actually, if we define 
the sequence of outcomes of the successive walks after 
a given number of steps as corresponding to a visit to 
node 2 or otherwise, we have a Bernoulli process with M 
observations and p = Pi(l, 2) to which we can associate 
the random variable Q(2) corresponding to number of 
times the agent reaches node 2. It can be shown that the 
random variable Q{2) follows the binomial distribution, 
with mean (Q(2)) = Mp and standard deviation crQ(2) — 
+ y^Mp{l — p). We can now define the mean frequency 
of accesses to node 2 as 



(i?(2)) = 



(g(2)) _ MFi(l,2) 



M 



M 



= Pi(l,2) = l/7 



which has standard deviation given as a-Q(2)/M (ob- 
serve that limM^oo <^q(2)I^ = 0). The mean period of 
accesses to node 2 is given as 



(r(2)} ^ 



M 



M 



(g(2)} MFi(l,2) Pi (1,2) 
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Now, let us consider the situation in which M self- 
avoiding walks start not only at node 1, but at all nodes 
(because the walks do not interact one another, they can 
be performed simultaneously or not). More specifically, 
M walks are initiated at each node of the network. We 
consider that each moving agent avoids only the nodes al- 
ready visited by itself, ignoring visits by the other nodes. 
Observe that several real systems of interest can be mod- 
eled by such walks, including the disease spreading and 
WWW surfing. After one step, all the moving agents 
starting from nodes 2 to 8 will necessarily be at node 1, 
so that the frequency of access to node 1 remains as be- 
fore and equal to 1/7. A completely different situation 
arises at the second step of the walks, in which case we 
have: 



FIG. 1; A simple, but fundamental, network (a)involving a 
central node and 7 surrounding nodes. Other simple networks 
considered for illustration of the accessibilities (b-d). Please 
refer to the text for discussion. 



r(M) 



0(2) 
M 



(1) 



After one step, the probability of reaching node 2 (or 
any of the other nodes 3 to 8), henceforth represented as 
Pi(l, 2), is formally defined as the following limit: 



P2(«,2) = l/6; 

where w G {2, 3, . . . , 8} and v G 3, 4, . . . , 8}. Because 
we have a total of N independent self-avoiding walks, it 
becomes reasonable to redefine the frequency of accesses 
to a node i, henceforth called inward accessibility of node 
i after h steps, as 



N 



q=l,q^i 



m)) 

M{N - 



N 

E 



q=l,q=il 



1) 



(2) 
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where the normaUzing factor (A'' — 1) is used instead of 
N because this corresponds to the maximum number of 
nodes which can reach node i after a single step (i.e. node 
i is naturally excluded from the normalization). There- 
fore, the inward accessibility of any node varies from 
to 1. Because of the non-linear dynamics of self-avoiding 
walks, the values of Ph{i,j) cannot be easily calculated 
in analytical fashion. In the current article, these prob- 
abilities are calculated computationally as described be- 
low. Observe also that = 1 a,nd, generally. 

The above definition of the inward accessibility mea- 
surement has as its main motivation the quantification 
of how accessible a given node i is by moving agents per- 
forming self-avoiding random walks, all of which starting 
from all nodes simultaneously. As the number of such 
walks increase (or more agents are liberated simultane- 
ously at each node), the inward accessibility of a node i 
will provide an indication of the intensity of visits to that 
node. 

Let us now address the motivation and definition of 
the outward accessibility of a node i after h steps of self- 
avoiding walks. Consider again the situation in Figure [U 
and let's assume that a self- avoiding random walk started 
at node 1. After the first step, any of the nodes 2 to 8 
will have an equal probability Pi(l,w) = 1/7 of being 
accessed. Such a non-preferential, equiprobable transi- 
tion from node 1 to the other nodes corresponds to the 
optimal situation favoring all the surrounding nodes to 
be accessed more uniformly and more speedily as succes- 
sive (or simultaneous) self-avoiding walks are performed 
starting at node 1. Such an effectiveness of accesses can 
be conveniently expressed [2l| in terms of the entropy 
(e.g. [l^) of the probabilities Pi(l,w). Therefore, in the 
above example we have this entropy Eh{i) is given as: 



Mi(l) = 



4/7 « 0.571 



N -I 

^-^^ = -ln{Q.2b)lln{7) « 0.712 

In order to avoid such a problem, in this work, we 
therefore define the outward accessibility as 



OAhii) 



exp{Eh{i)) 
N - 1 



(3) 



Observe that a node having just one immediate neigh- 
bor will have the smallest entropy value OA = 1/{N — 1) 
(we are not considering networks with a single isolated 
node). 

Now, going back to the example in Figure [ijb) , we 
have: 



Ei{l) = -;n(0.25) 



so that 



OA.(l) = ""(-Mf-^'S)) , 4/7 

i.e. the inward and outward accessibilities will be equal 
for the specific situation involving equiprobability transi- 
tions from a central node. Also, we necessarily have that 
< I A < 1. A situation involving non-equiprobable 
transitions is shown in Figure [TJc). In this case: 



JV 



i?i(l) = -^Pi(l,<7)Zn(Pi(l,(7)) 



Though such a feature provides a sound basis for quan- 
tifying the diversity of the walks originating at i, explored 
partially in [2lj , it has the shortcoming of exhibiting val- 
ues comprehended between and ln{N — 1), and not be- 
tween and 1 as is the case with the above defined inward 
accessibility. The possibility of normalizing Eh{i) by di- 
viding it by ln{N — 1) presents the inconvenient that the 
inward and outward accessibilities involving transitions 
with the same probabilities would have different values. 
In order to illustrate such a problem, refer to the situa- 
tion depicted in Figure [Tfb), where the central node 1 is 
connected to 4 surrounding nodes (1, 2, 3 and 4), while 
other nodes are connected to node 4. The inward acces- 
sibility and diversity entropy of node 1 can be calculated 
as 



IA,{1) 



N - 1 



O^i(l) 



ea;p(0.1/n(0.1) -f 0.9;n(0.9))(2) 



2/2 = 1 
« 0.692 



Such values properly reflect the fact that node 1 in 
the network in Figure [ijc) has greater (actually optimal) 
potential to receive accesses from other nodes after one 
step than to access other nodes after that same number 
of steps. Actually, in the limit when Pi (1,2) becomes 
very small, we have 



lim 

(1^0 exp{aln{a) + (1 — a}ln{l — a)) 



= 0.5 



(4) 



which, as could be expected, corresponds to the situa- 
tion illustrated in Figure [TJd). 

It is interesting to observe that the quantity 
\/exp{—Eh{i)) can be understood, from the perspec- 
tive of outward accessibility, as the 'equivalent number 
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of nodes' attached to the reference node i, a real value 
instead of an integer. In other words, the outward ac- 
cessibility of node 1 in the network in Figure [T] would 
be equivalent to the accessibility of that node connected 
to l/exp{—Ei(l)) 'fractionary nodes' with equiprobable 
transitions exp{~Ei{l)). In that case, we have 1.384 
equivalent nodes and transition probabilities 0.722. 

It should be observed that the above definitions of the 
inward and outward accessibilities are independent of the 
type of walk under consideration and can be easy and 
immediately extended to all other types of random walks 
and Markovian systems, linear or not, provided the tran- 
sition probabilities P/j(i,j) can be estimated or calcu- 
lated. Actually, such probabilities can be understood as 
corresponding to the strengths of respectively defined vir- 
tual networks, in the sense that each probability Ph{i,j) 
can be interpreted as a virtual link [3(| between nodes i 
and j at distance h. Consequently, it becomes reasonable 
to speak of nodes with high probabilities as hubs for the 
respectie h. 

Because of the non-linear dynamics of the self-avoiding 
walks, the transition probabilities Ph{i,j) between all 
pairs of distinct nodes are difficult to be analytically cal- 
culated. In this work, we adopt the simple computational 
procedure described in [2l| . which involves performing 
several self-avoiding walks from each node and accumu- 
lating the relative frequencies of access to each visited 
node in terms of the number of steps h. 



III. RESULTS AND DISCUSSION 



A representative network has been chosen for each of 
the 6 considered network types in order to illustrate the 
estimation and interpretation of the inward and outward 
accessibilities of nodes. All networks had (k) w 6 and 
N — 100, except for the GG case, which has = 91 
nodes. Because the length of the self-avoiding walks 
were observed to be almost always equal to the maximum 
number of steps considered in this work (i.e. 10 steps), 
the obtained results are reasonably devoid of border ef- 
fects (i.e. most of the walks did not reach the extremity 
nodes). This does not mean, however, that the following 
results will not exhibit specific scaling effects with the 
network size A'^. 

Figures [J] to [5] show the respective measurements 
for the ER, BA, WS, GG, PN and PA networks. The 
y— axes corresponds to the number of steps along the 
self-avoiding walks (i.e. h). The gray intensity in these 
images indicate higher values of the measurements, nor- 
malized with respect to the maximum and minimum val- 
ues for the sake of better visualization. A total of 2000 
self-avoiding random walks were performed for each node 
by the algorithm for numeric estimation of the transition 
probabilities Ph{i,j)- 

The first important result is that the inward accessi- 
bility values were verified to be much smaller than the 



respective outward accessibilities, as can be verified from 
Table m Observe that the mean I A values are necessarily 
0.001 — l/N because of the normalization of the proba- 
bilities. Part of the reasons for the asymmetry of accesses 
shown in this table is readily clear from the example in 
Figure [TJa): after two steps along the walks, any of the 
surrounding nodes will have OA2 = 6/7 and IA2 = 1/7. 
In this case, the marked asymmetry is caused by the dif- 
ference of degrees between the central and the surround- 
ing nodes as well as the fact that the surrounding nodes 
are attached to the central node only. 

Because the regularity of networks is always a char- 
acteristic affecting several of the topological properties, 
let us consider how the accessibilities behave in a lattice 
network such as that shown in Figure [D Such a network 
is highly regular with respect to the node degrees, except 
at its borders. Let us study the central node. At step 
h = l,we have OAi = 4/48 and lAi = 1/48. For h = 2, 
we have OA2 = 8/48 and IA2 = 1/48. Observe that the 
asymmetry tends to increase strongly with h. Similar 
situations are verified for the other nodes in the network. 
Note that in the case of the network in Figure [21 the ac- 
cessibilities to the central node for h — 1 and h = 2 are 
not affected by border effects. 

A little reasoning reveals that the marked asymme- 
try between the outward and inward accessibilities in 
many networks is ultimately a consequence of the de- 
gree regularity and high average degree, implying that 
usually each node in a network tends to be connected to 
many additional distinct nodes as h increases. Observe 
also a striking difference between the situations defined 
by Pi{central, surround) in the networks in Figure a) 
andO while the surrounding nodes in the former case can 
only move to the central node after one step (implying the 
central node to have maximum inward accessibility), in 
the latter situation the nodes surrounding the reference 
node can move to many other places (therefore decreas- 
ing substantially the inward accessibility of the central 
node). This effect, which acts as a 'diode' on the flow of 
accesses, tends to increase substantially with h and is the 
main cause of accessibility asymmetry in typical complex 
networks. 

Additional insights about the nature of asymmetries 
in specific cases can be gathered by inspecting Figure [31 
which show the PGA outward (a) and inward (b) acces- 
sibilities for the considered GG network. These measure- 
ments were obtained by applying principal component 
analysis (PGA) [H, [3l[ considering all number of steps, 
so that the resulting principal variables corresponded 
to the most relevant uncorrelated linear combination of 
measurements. First, observe that a group of nodes (in 
green) at the center of the graph resulted with the largest 
OA values (Figure [31Ia)). This is because such nodes 
have more neighbors along the initial steps of the self- 
avoiding walks than the nodes which resulted in smaller 
OA values, which correspond to the 'border' of the graph. 
Observe that the OA values can therefore be used to de- 
fine the border of the network as corresponding to the 
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Network 


OA 


lA 




mean ± st. dev. 


min 


max 


mean ± st. dev. 


min 


max 


ER 


0.73 ± 0.30 


0.001 


0.73 


0.010 ±0.003 


0.001 


0.027 


BA 


0.64 ±0.26 


0.02 


0.26 


0.010 ±0.008 


0.001 


0.096 


WS 


0.52 ±0.28 


0.03 


0.94 


0.010 ±0.001 


0.004 


0.016 


GG 


0.28 ±0.17 


0.001 


0.76 


0.010 ± 0.004 


0.002 


0.035 


PN 


0.77 ±0.32 


0.50 


0.98 


0.010 ±0.001 


0.008 


0.011 


PA 


0.71 ±0.31 


0.20 


0.94 


0.010 ±0.003 


0.003 


0.023 



TABLE I: The mean, standard deviations, minimum and maximum values of the outward and inward accessibilities obtained 
for each of the networks, considering all the 10 steps. Because of the normalization of the transition probabilities, the mean of 
I A is always equal to 0.01. 
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FIG. 2: A lattice network with A'^ = 49 nodes. The studied 
node (central) is marked with an asterisk, and the nodes which 
can be visited by the central node after two steps are shown 
in black. 



nodes with the smallest OA values. Now, let us anal- 
yse the I A values shown in Figure [Hb) . An opposite 
tendency is clearly observed in this case, with the 'bor- 
der' nodes yielding the highest I A values. This is be- 
cause such groups of nodes tend to act as traps for the 
self-avoiding walks (recall that the moving agents remain 
at the termination nodes). Observe also that the nodes 
leading to such traps (e.g. 3, 75 and 82) tend to have 
the smallest values of lA, as they do not accumulate 
the agents. The presence of traps, i.e. nodes or group of 
nodes loosely connected with the rest of the network, con- 
tributes greatly to decreasing the I A values at the other 
nodes. In the specific case of the type of GG network 
considered in this work, such border (in the topological 
sense) nodes and groups of nodes tend to result at the 
spatial border of the graph. 

A series of interesting results and insights can be in- 
ferred from Figures 2] to [51 which shows the evolution 
(bottom-up, along the vertical direction in the figures) 
of the OA, I A and outward/inward ratios along the 10 
steps (y— axes) for each of the network nodes (along the 
a;— axes). We discuss such results from the general to 



the specific perspective as follows. To start with, OA in- 
creases with h in all cases. In agreement with previous 
results considering the diversity entropy signatures [2l| . 
the increase of OA with h tended to be more gradual for 
the WS and GG networks, possibly for the same reasons 
identified in 12 1| . namely the fact that pairs of nodes 



in these two types of networks tend to be less adjacent 
regarding longer paths fs^ . I A tended to remain nearly 
constant with h for ER, WS, PN and PA. This is a sur- 
prising result, because it indicates that the inward ac- 
cessibility does not seem to depend on the length of the 
walks. As we will discuss later in this section, the I A val- 
ues are strongly positively correlated with the individual 
node degree in most cases, especially for /i > 3. In the 
case of GG (see Figure (TJb)), the I A values tended to 
increase, decrease or to present a peak along the steps h. 
Such a marked difference between this model and the 5 
other types of networks suggests that the invariance of I A 
with h in the latter models seems to be related to some 
intrinsic feature of the GG type of network, possibly its 
strictly local node adjacency and lack of bypasses short- 
ening the distance between nodes, so that the I A values 
in this case are more dependent on the local connectivity 
structure of the network. 

As shown in FigurelHl the I A in the case of the BA net- 
work was dominated by the hub corresponding to node 2, 
which has degree 38. Interestingly, the I A values tended 
to oscillate with h for this node, decreasing for the largest 
values of h. Another interesting result regards the fact 
that OA tended to get more uniform among the nodes for 
the highest values of h, except for the GG model. This 
is a consequence of the fact that the OA values tend to 
increase steeply for most models (except WS and GG) 
towards plateaux of similar heights for most nodes in the 
networks. Therefore, after 2 or 3 time steps, most nodes 
in the ER, BA, WS, PN and PA will tend to access a 
large portion of the network nodes. The long range con- 
nections and path-adjacencies in these networks seem to 
play an important role in this effect. In the case of the 
GG structure, some nodes (e.g. nodes 4, 6, 8, 9, etc.) 
resulted with markedly higher values of OA along the 
highest values oih- see Figure [71 This implies that some 
nodes in the network acquire an especially important role 
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in accessing nodes as the number of steps increase in this 
type of networks. 

Unhke OA, the values of I A tend to vary little with h 
for each node, except for the WS and PN cases, in the 
sense that substantial similarity is observed among the 
I A values of several nodes. As a matter of fact, I A be- 
comes intensely uniform for the largest values of h for the 
PN network (see Figure [8]) , corroborating the enhanced 
uniformity of this type of network (see also [21,, i22, ;2^). 
The few nodes in Figure [5] presenting smaller values of 
I A have been verified to correspond to the initial and 
final nodes of the self-avoiding walks which are used to 
construct that type of network [221 . 

Because the inward accessibilities are largely smaller 
than the respective outward counterparts, the values of 
the ratios outward/inward tended to mostly reflect the 
latter measurement, except mainly for the PN case. More 
specifically, the outward/inward ratios in this network 
indicates that the extremity nodes used to build this type 
of network tend to present a tendency to access more than 
being accessed, a trend which increases with h. 

Having investigated the overall properties of the in- 
ward and outward accessibilities for examples of the 6 
considered models, we now turn our attention to the par- 
ticularly relevant problem of trying to predict the dynam- 
ical property of accessibility by considering exclusively 
topological features of the networks. Observe that this 
issue is intrinsically related to the structure-dynamics 
paradigm in complex networks research. Though such 



an analysis can be performed while considering several 
network measurements [23^, in the current article we fo- 
cus our attention on correlations between the individual 
node accessibility and the respective degree k. A similar 
investigation targeting the correlation between linear dif- 
fusion (simulated by traditional random walks) and node 
degree has been reported in . That work showed that 
though full positive correlation is verified between visits 
received by nodes and the respective degree in the case 
of undirected, connected networks, such a correlation is 
generally broken in the case of directed networks. A suf- 
ficient condition, namely that the in and outdegree of all 
nodes be identical, was identified in that work. The con- 
sideration of non-linear dynamics implemented by self- 
avoiding networks in the current article naturally breaks 
the full correlation between node accessibility and degree. 
Yet, a strong positive correlation between the inward ac- 
cessibility and node degree has been verified for several 
of the 6 networks studied in the current work, especially 
for h larger than 3 or 4 steps. 

Figures [TU] to [H] illustrate the scatterplots of OA x k 
and lAxk respectively to the ER, BA, WS, GG, PN and 
PA networks and h — 1 (the initial step), 5 (the middle 
step) and 10 (the final steps in our self-avoiding walks). 
Figure [1^] shows the evolution of the Pearson correlation 
coefficients for the correlations I A x k and OA xkin terms 
of h for all networks. Recall that the Pearson correlation 
coefficient varies from -1 (extreme negative correlation) 
to 1 (extreme most positive correlation), with null value 
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FIG. 4: The outward accessibility (a), inward accessibility (b) and the outward/inward ratio (c) for the 10 initial steps obtained 
for the considered ER network with TV — 100 nodes and (fc) = 6. 
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FIG. 5: The outward accessibility (a), inward accessibility (b) and the outward/inward ratio (c) for the 10 initial steps obtained 
for the considered BA network with A'^ — 100 nodes and (fc) — 6. 



indicate that the two measurement are uncorrelated. 

The first important conclusion regarding the relation- 
ship between the outward accessibility and the individual 
node degree concerns the fact that a full positive corre- 
lation is obtained between these two features for h = 1 



(see Figures m to Eland Figure fTB]) . However, such a cor- 
relation decreases steeply with h, and becomes almost 
irrelevant after 4 or 5 steps (see Figure [16]) . These re- 
sults imply that the outward accessibility of each node 
can only be used to predict the respective outward ac- 
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FIG. 6: The outward accessibility (a), inward accessibility (b) and the outward/inward ratio (c) for the 10 initial steps obtained 
for the considered WS network with N = 100 nodes and (k) — 6. 
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FIG. 7: The outward accessibility (a), inward accessibility (b) and the outward/inward ratio (c) for the 10 initial steps obtained 
for the considered GG network with A'^ = 91 nodes and (k) — 6. 



cessibility for the 3 or 4 initial steps. 

Interestingly, an opposite trend has been observed for 
the correlation between the inward accessibility and the 
degree, which is generally high at h — I and tends to 
increase further with h for most cases, except WS and 



GG (see Figure [TB|) . This result implies that the inward 
accessibility of a node can be predicted with remarkable 
accuracy from the respective node degree, especially for 
relative high values of h. This is surprising because the 
higher the value of h, the more global the dynamics is un- 
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FIG. 8: The outward accessibility (a), inward accessibility (b) and the outward/inward ratio (c) for the 10 initial steps obtained 
for the considered PN network with N = 100 nodes and (k) = 6. 
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FIG. 9: The outward accessibility (a), inward accessibility (b) and the outward/inward ratio (c) for the 10 initial steps obtained 
for the considered PA network with N = 100 nodes and (k) = 6. 



folded along the topology of the network. The surprise 
therefore resides in the fact that the node degree, an in- 
trinsically local property of the network connectivity, can 
be used to predict accurately (except for WS and GG) 
the global transient dynamics of the self-avoiding walks 



for most values of h. Such a result has particularly im- 
portant practical implications. For instance, in a system 
properly modeled by the self-avoiding walks unfolding on 
complex networks similar to ER, BA, PN and PA, it is 
possible to predict the number of visits to each node in 
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FIG. 10: The scatterplots of the outward and inward accessibihties agains the node degree in the ER network, with respective 
Pearson correlation coefficients. 
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FIG. 11: The scatterplots of the outward and inward accessibilities against the node degree in the BA network, with respective 
Pearson correlation coefficients. 
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FIG. 12: The scatterplots of the outward and inward accessibihties against the node degree in the WS network, with respective 
Pearson correlation coefficients. 



(a) (OA X k) at step 1 for GG, PeaiHjn= 1 



0.14 
0.12' 
0.10 
0.08 
0.06 
0.04 
0.02 




10 



12 
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FIG. 13: The scatterplots of the outward and inward accessibilities against the node degree in the GG network, with respective 
Pearson correlation coefficients. 



13 



(a) (OA X k) at step I for PA, Pearson= 1 



(b) (OA X k) at step 5 for PA, Pearson= .55 



(c) (OAXk) at step 10 for PA, Pearson= .17 




945e-3 
94Ie-3 
9ile-} 
933e-3 
929e-3 
925e-3 
92Ie-3 
917e-3 
913e-3 
909e-3 
905e-3 







(f) ([A X fc) at step 10 for PA, Pearscm= 1 



20e-3 
18e-3 
16e-3 
14e-3 
12c-3 
lOe-3 
8e-3 
6e-3 
4e-3 



10 



12 



FIG. 14: The scatterplots of the outward and inward accessibihties agains the node degree in the PA network, with respective 
Pearsotn correlation coefficients. 
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FIG. 15: The scatterplots of the outward and inward accessibilities against the node degree in the PN network, with respective 
Pearson correlation coefficients. 
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terms of its degree. 

The above result has major impHcations for several 
real-world problems. For instance, in the case of disease 
spreading or node attack, the more susceptible nodes can 
be easily identified and protected. If the dynamics is used 
to model neuronal or dissemination of cortical activation, 
the more active nodes can be immediately identified from 
their degrees. The strong correlation between the inward 
accessibility and the individual node degree also paves 
the way for interfering in the local network topology (e.g. 
changing the degree of specifically critical nodes) in or- 
der to address specific needs. For instance, in the case of 
WWW surfing, a site (i.e. a node) can get more accessed 
provided it can motivate additional links (however, note 
that this type of network would typically be directed, so 
that such an application would require additional studies 
on self-avoiding walks on that type of asymmetric connec- 
tivity) . The correlation between inward accessibility and 
degree verified for most of the networks is strongly related 
to the fact that, once the walks have become spread along 
the network as h increases (the ER, BA, FN and PA net- 
works tend to have strong and regular path-adjacency), 
most edges tend to exhibit similar rates of accesses, so 
that the final inward accessibility is ultimately defined 
by the individual node degree, which taps visits with an 
intensity which is almost linear with the number of in- 
coming edges. The strategy of enhancing the I A corre- 
lation by adding new connections with other nodes will 
not generally work in the GG case because in that case 
a node typically receives connections only from nearby 
nodes (regarding spatial distance) , so that it will tap ac- 
cesses from nodes with similar — and possibly low — 
accessibility. Ultimately, it is such a possibility of segre- 
gation of levels of accessibility which contributes to the 
break of the I A x k correlation in the case of the GG 
(and also to a certain extent in the WS) model. 

Interestingly, the importance of uniformity of edge ac- 
tivation for the inward correlation also explains why it 
is so difficult to predict the outward accessibility from 
the degree larger values of h: the problem is that the 
dissemination of the accesses tends to become uniform 
as the walks unfold along longer lengths irrespectively of 
the local connectivity along the initial steps of the walks. 
These results and interpretations suggest investigations 
of the accessibility to edges, in addition to the studies 
targeting nodes reported in the current work. 

Still considering the positive correlation between the 
inward accessibility and degree. It can be verified from 
FigureOthat the correlation in the scatterplot for h = 10 
suggests a saturation of that relationship, manifested by 
the sigmoid bending at the right-hand side of the curve. 
Another interesting issue regards the reasons for the lack 
of correlation between the inward accessibility and the 
degree observed for the GG model. Strictly speaking, 
some significant correlation is observed for /i = 1 to 3 
or 4 (see Figure [TH]) , but it tends to fade away for larger 
values of h. One possible explanation is that the more 
local adjacency between nodes in this type of network (a 



consequence of spatial adjacency, which tends to present 
communities for the considered size — see Figure ^ , im- 
plies that the flow of accesses along the edges of the net- 
work for large values of h will not be so uniform, making 
it impossible to predict the inward accessibility in terms 
of the individual node degrees. 

If the GG network implied the worst predictability, the 
PA allowed the opposite feature, exhibiting Pearson cor- 
rellation coefficient for I A x k very close to 1 except for 
h = 1. This result seems to be related to the fact that 
paths of several lengths are guaranteed to exist among 
several of the nodes in this type of networks as a conse- 
quence of the way in which they are generated, i.e. by 
transforming the star connectivity into path connectivity. 



IV. CONCLUDING REMARKS 

One particularly interesting means to investigate the 
relationship between the dynamics and structure in com- 
plex networks involves the study of non-linear random 
walks such as the self-avoiding walks. Unlike traditional 
walks, that type of walk is more purposive and effective 
in the sense of economically visiting nodes, i.e. with the 
minimal possible number of edges. Because such a kind 
of walk cannot repeat nodes or edges, it implies in a more 
effective diffusion of visits away from the original node, 
without the returns implied by the traditional random 
walks. In addition, unlike their traditional counterparts, 
self-avoiding walks are necessarily finite. Despite such 
important and useful features, self-avoiding walks have 
received relatives little attention from the complex net- 
works literature. Even less frequently studied is the tran- 
sient dynamics in such systems. 

Another important aspect regarding the study of dy- 
namics in complex networks regards the accessibility of 
a node. This property is critical because it underlies a 
series of important theoretical and real- world problems, 
including disease spreading, WWW navigation, protein 
interaction, transportation, urban planning, communica- 
tions, distributed computing, cortical network activity, 
neuronal networks (where the refractory period of neu- 
ronal firing can be modeled by the self-avoiding crite- 
rion), amongst many other cases. 

The current article has brought together these two rel- 
evant issues, namely self-avoiding walks and node in ward 
and outward accessibility. More specifically, we build up 
from previous results [211 in order to propose two sound 
and meaningful definitions of the accessibilities. There- 
fore, by outward accessibility we quantify the potential 
of each node in a given network to be accessed by self- 
avoiding walks originating at all other network nodes. 
The inward accessibility has an analogue meaning re- 
garding the accesses received by each specific node. Both 
such measurements can be calculated from the transition 
probabilities between any pair of distinct nodes for each 
considered number of steps. In this work, these proba- 
bilities have been estimated by using a simple algorithm 



15 




FIG. 16: The Pearson correlation coefficients obtained between the outward and inward accessibihties and the node degrees 
for each of the considered networks. 



which consists exactly in implementing the considered 
dynamics, i.e. performing a large number of self-avoiding 
walks from each node. The particularly relevant issue of 
trying to predict the dynamical property of accessibility 
by considering exclusively topological features of the net- 
works (in the case of the present work the individual node 
degree) has also been addressed. This issue is at the very 
center of the structure/dynamics relationship paradigm. 

Because of the quantity of results, the reader is kindly 
requested to refer to Section Hill for a complete listing and 
discussion of the results. Only the most relevant results 
and interpretations are listed in the following. The first 
result which stands out regards the fact that the outward 
accessibility is invariably much larger than the inward ac- 
cessibility for most nodes. This has been verified to be a 
consequence of the relatively large degrees and regularity 



observed for most of the considered networks, which en- 
hances the possibilities that the outgoing walks diverge 
among the network. As a matter of fact, observe that the 
outward accessibility presents an analogy with the Lya- 
punov coefficient in dynamics systems, in the sense that 
it quantifies the divergence of the dynamics along the 
consecutive steps of the walks. Another particularly rel- 
evant result obtained in this work concerns the fact that 
most accessibility values for each individual node tend to 
increase with h. A series of results and insights were 
additionally obtained regarding the application of the 
measurements to 6 distinct types of complex networks, 
namely Erdos-Renyi, Barabasi- Albert, Watts- Strogatz, a 
simple geographical type of network, as well as two knit- 
ted networks (PN and PA). Out of these models, the GG 
case led to the most different dynamics, as a consequence 
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of its intense local adjacency. The accessibilities were 
found to increase more abruptly along h for the ER, BA, 
PN and PA structures. Except for the GG model, the 
outward accessibility was moderately correlated with the 
node degree along the initial steps for all models, fading 
away after 3 or 4 steps. Surprisingly, the inward acces- 
sibility tended to be high and to increase further with 
the number of steps (except of the WS and GG case), 
implying that such a global dynamical property can be 
accurately predicted from the degree of individual nodes. 
This important result is related to the uniformization of 
the flow accesses among the edges typically achieved for 
larger values of h. 

Several are the possibilities for future investigations 
motivated by the currently reported concepts, methods 
and results. To begin with, it would be interesting to con- 
sider other values of average degree and network sizes, 
though this will require considerable computational re- 
sources. Because of the intrinsic importance of accessi- 
bility, it would be particularly interesting to characterize 
the inward and outward accessibility of real-world net- 
works, such as those related to disease spreading, trans- 
portation, biological interactions, amongst many others. 
Also, as the definition of accessibilities proposed in this 
article are independent of the type of network and ran- 
dom walk, it would be easy and immediate to extend 
such analyses to additional theoretical network models 
and less conventional random walks such as those de- 
scribed in T^, including direct and weighted structures. 
Particularly promising is the application of random walks 
formulation to model and study synchronization in linear 
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